Evolutionary history of the mtDNA 9-bp deletion in Chinese populations and its relevance to the peopling of east and southeast Asia

In total, 1218 Chinese from twelve ethnic groups and nine Han geographic groups were screened for the mtDNA 9-bp deletion motif. The frequency of the 9-bp deletion in all samples was 14.7% but ranged from 0% to 32% in the various ethnic groups. Three individuals had a triplication of the 9-bp segment. Phylogenetic and demographic analyses of the mtDNA hypervariable segment 1 (HVS1) sequences suggest that the 9-bp deletion occurred more than once in China. The majority of the Chinese deletion:haplotypes (about 90%) have a common origin as a mutational event following an initial expansion of modem humans in eastern Asia. Other deletion haplotypes and the three haplotypes with a 9-bp triplication may have arisen independently in the Chinese, presumably by replication error. HVS1 haplotype analysis suggests two possible migration routes of the 9-bp deletion in east and southeast Asia. Both migrations originated in China with one route leading to the Pacific Islands via Taiwan, the other to southeast Asia and possibly the Nicobar Islands. Along both routes of peopling, a decrease in HVS1 diversity of the mtDNA haplotypes is observed. The "Polynesian motif (16217T/C, 16247A/G, and 16261C/T)" and the 16140T/C, 16266C/A, or C/G polymorphisms appear specific to each migration route.

Reconstructing the evolutionary history of China: A caveat about inferences drawn from ancient DNA

The decipherment of the meager information provided by short fragments of ancient mitochondrial DNA (mtDNA) is notoriously difficult but is regarded as a most promising way toward reconstructing the past from the genetic perspective. By haplogroup-specific hypervariable segment (HVS) motif search and matching or near-matching with available modem data sets, most of the ancient mtDNAs can be tentatively assigned to haplogroups, which are often subcontinent specific. Further typing for mtDNA haplogroup-diagnostic coding region polymorphisms, however, is indispensable for establishing the geographic/genetic affinities of ancient samples with less ambiguity. In the present study, we sequenced a fragment (similar to 982 bp) of the mtDNA control region in 76 Han individuals from Taian, Shandong, China, and we combined these data with previously reported samples from Zibo and Qingdao, Shandong. The reanalysis of two previously published ancient mtDNA population data sets from Linzi (same province) then indicates that the ancient populations had features in common with the modem populations from south China rather than any specific affinity to the European mtDNA pool. Our results highlight that ancient mtDNA data obtained under different sampling schemes and subject to potential contamination can easily create the impression of drastic spatiotemporal changes in the genetic structure of a regional population during the past few thousand years if inappropriate methods of data analysis are employed.

The ADH1B Arg47His polymorphism in East Asian populations and expansion of rice domestication in history

Background: The emergence of agriculture about 10,000 years ago marks a dramatic change in human evolutionary history. The diet shift in agriculture societies might have a great impact on the genetic makeup of Neolithic human populations. The regionally restricted enrichment of the class I alcohol dehydrogenase sequence polymorphism (ADH1BArg47His) in southern China and the adjacent areas suggests Darwinian positive selection on this genetic locus during Neolithic time though the driving force is yet to be disclosed. Results: We studied a total of 38 populations (2,275 individuals) including Han Chinese, Tibetan and other ethnic populations across China. The geographic distribution of the ADH1B*47His allele in these populations indicates a clear east-to-west cline, and it is dominant in south-eastern populations but rare in Tibetan populations. The molecular dating suggests that the emergence of the ADH1B*47His allele occurred about 10,000 similar to 7,000 years ago. Conclusion: We present genetic evidence of selection on the ADH1BArg47His polymorphism caused by the emergence and expansion of rice domestication in East Asia. The geographic distribution of the ADH1B*47His allele in East Asia is consistent with the unearthed culture relic sites of rice domestication in China. The estimated origin time of ADH1B*47His allele in those populations coincides with the time of origin and expansion of Neolithic agriculture in southern China.

银杉的分子谱系地理学研究

银杉（Cathaya argyrophylla）是中国特有的濒危裸子植物，孤立分布于我国亚热带山地。虽然以往等位酶和RAPD标记的研究表明，银杉群体的遗传多样性水平很低而群体间的遗传分化很高，但迄今对该种群体的动态变化以及物种的进化历史仍缺乏了解，包括影响群体遗传结构的因素以及物种可能的避难所等，对如何有效地保护和恢复银杉群体仍缺乏科学依据。本文根据8个核基因片段和2个线粒体片段的序列数据，运用群体遗传学和谱系地理学方法，探讨了银杉在DNA水平上的多样性和群体的动态历史，探讨了影响银杉特殊的群体遗传结构的各种因素，并推测了银杉第四纪冰期的避难所，对银杉花粉活力及其变异进行检测。在此基础上，提出了保护和恢复银杉群体的具体策略和措施。主要研究结果如下： 　　1. 核甘酸多态性和群体遗传结构 　　从101个核基因片段中筛选了8个适用于银杉的片段，对来自四个地区15个群体共86个个体的胚乳总DNA进行了扩增和测序。8个核基因座位的平均核苷酸多态性（θs＝0.0022，πs＝0.0027）显著低于其它松杉植物遗传多态性的多座位估计值。四个地区中，大瑶山（DY）的多态性最高（θs＝0.0026，πs＝0.0027），八面山（BM）最低（θs＝0.0014，πs＝0.0018），大娄山（DL）和越城岭（YC）介于二者之间。大多数座位内的多态位点间紧密连锁，座位间的连锁只在八面山地区检测到。AMOVA分析表明显著性的多态性比例存在于地区间（20.05％）和地区内群体间（9.37％）。FST分析也显示群体间（FST＝0.294）和地区间（FST =0.131-0.319）存在显著的遗传分化。推测伴随着瓶颈效应而出现的遗传漂变及其有限的基因流（Nm=1.2）是导致银杉群体低水平多态性和群体间强烈分化的主要原因。 　　2．谱系地理学分析 　　利用2个线粒体片段（nad1和nad4）序列以及高变异量的核2009片段序列对银杉的谱系地理进行了探讨。2个线粒体片段的多态位点组合成3种单倍型，将银杉分成大娄山（DL）、八面山（BM）以及越城岭（YC）和大瑶山（DY）3个地区，地区间的单倍型完全不同（GST=1.0），结合核基因呈现的群体遗传结构，推测现存银杉群体至少来源于4个冰期避难所，相当于银杉现存的4个相互隔离的地区。2009座位上12个变异位点组合成8种单倍型，位于单倍型谱系内部节点的4种祖先单倍型分布广泛、出现频率最高，其它7个核基因座位具有类似谱系结构。遗传距离和地理距离没有相关性，NST (0.138)与GST (0.134)没有显著性差异，说明现存的银杉群体是相对较近的时间内片断化的产物。2009片段分离位点的失配分布（mismatch distribution）呈双峰和多峰，表明银杉群体没有经历近期的扩张，与古生物学研究证据相吻合。 　　3. 银杉的基因杂合性和花粉生命力 　　利用2009和cad两个核基因片段，采用多胚乳序列法得到总的杂合体比率为64％（2009）和60％（cad），说明银杉群体中存在高比例的杂合体。大娄山地区的杂合体比率是八面山地区的2倍。银杉杂合体比率的高低可能与其遗传多态性有关，也可能是自然选择的结果。采用TTC染色法对银杉的花粉生命力进行了测定，在干燥低温条件下银杉花粉的活力很稳定，保存一年后有活力的花粉数仍高达80％以上。通过对来自两个地区（DY和YC）7个群体共16个银杉个体花粉活力的测定发现，银杉有活力花粉比例高达93.3％，与其它裸子植物相当。花粉活力在地区间和群体间存在显著差异，花坪地区（95.2%）的花粉活力高于大瑶山（91.3%）。花粉活力在群体内个体间差异不显著。 　　4．银杉的进化历史及其保护 　　银杉的低水平的遗传多样性和独特的群体遗传结构对我们推测其冰川期避难所提供了重要依据。本研究在银杉4个孤立分布区发现了彼此不重叠的线粒体单倍型，同时核基因表现出了4个地理群，说明随着第四纪冰期气候的波动和银杉分布范围的片段化，原来广布的群体逐步萎缩，最后被保留在位于西部大娄山（DL），东部八面山(BM)，中南部越城岭(YC)和南部大瑶山(DY) 4个相互隔离的避难所。银杉独特的群体结构和动态历史对进一步制定相应的保护措施具有重要参考价值。由于遗传多态性很低，群体又小，几乎所有现存的银杉群体都面临由于随机事件导致的物种灭绝。更严重的是，当前该物种的适生环境不断恶化和片段化，以及异常低的生殖和存活率导致银杉与其它物种竞争能力很低。因此，除目前采取的原地保护策略外，迁地保护应给予优先考虑。此外，采用传统的控制授粉策略（在遗传上有显著差异的群体间开展人工杂交）是丰富其遗传多态性、恢复衰退群体的可行措施之一。

Molecular phylogeny of wood mice (Apodemus, Muridae) in East Asia

Sequences of the mitochondrial cytochrome b (1140 bp) and nuclear IRBP (1152 bp) genes were used to assess the evolutionary history of Apodemus, using the complete set of Asian species. Our results indicate that speciation in Asia involved three radiations, which supports an earlier study. The initial radiation yielded A. argenteus (Japanese endemic), A. gurkha (Nepalese endemic), and the ancestral lineage of the remaining Asian species. This lineage subsequently diverged into four groups: agrarius-chevrieri (agrarius group), draco-latronum-semotus (draco group), A. peninsulae, and A. speciosus (Japanese endemic). The final step consisted of divergence within two species groups as a consequence of the geography of the Yunnan-Guizhou plateau and Taiwan. The ecological ability of two Apodemus-species to inhabit one locality via niche partitioning likely drove the second radiation and shaped the basic geographical pattern seen today: A. argenteus and A. speciosus in Japan, A. agrarius and A. peninsulae in northern China, and the A. agrarius and A. draco groups in southern China. The three radiations are estimated to have occurred 7.5, 6.6, and 1.8-0.8 Mya respectively, using the IRBP clock, based on rat-mouse divergence 12 Mya. (C) 2003 The Linnean Society of London.

Molecular phylogeny of Sinocyclocheilus (Cypriniformes : Cyprinidae) inferred from mitochondrial DNA sequences

More than 10 species within the freshwater fish genus Sinoncyclocbeilus adapt to caves and show different degrees of degeneration of eyes and pigmentation. Therefore, this genus can be useful for studying evolutionary developmental mechanisms, role of natural selection and adaptation in cave animals. To better understand these processes, it is indispensable to have background knowledge about phylogenetic relationships of surface and cave species within this genus. To investigate phylogenetic relationships among species within this genus, we determined nucleotide sequences of complete mitochondrial cytochrome b gene (1140 bp) and partial ND4 gene (1032 bp) of 31 recognized ingroup species and one outgroup species Barbodes laticeps. Phylogenetic trees were reconstructed using maximum parsimony. Bayesian, and maximum likelihood analyses. Our phylogenetic results showed that all species except for two surface species S. jii and S. macrolepis clustered as five major monophyletic clades (I, II, III, IV, and V) with strong supports. S. jii was the most basal species in all analyses, but the position of S. macrolepis was not resolved. The cave species were polyphyletic and occurred in these five major clades. Our results indicate that adaptation to cave environments has occurred multiple times during the evolutionary history of Sinocyclocheilus. The branching orders among the clades I, II, III, and IV were not resolved, and this might be due to early rapid radiation in Sinocyclocheilus. All species distributed in Yunnan except for S. rhinocerous and S. hyalinus formed a strongly supported monophyletic group (clade V), probably reflecting their common origins. This result suggested that the diversification of Sinocyclocheilus in Yunnan may correlate with the uplifting of Yunnan Plateau. © 2005 Published by Elsevier Inc.

Molecular evolution of growth hormone gene family in old world monkeys and hominoids

Growth hormone is a classic molecule in the study of the molecular clock hypothesis as it exhibits a relatively constant rate of evolution in most mammalian orders except primates and artiodactyls, where dramatically enhanced rate of evolution (25-50-fold) has been reported. The rapid evolution of primate growth hormone occurred after the divergence of tarsiers and simians, but before the separation of old world monkeys (OWM) from new world monkeys (NWM). Interestingly, this event of rapid sequence evolution coincided with multiple duplications of the growth hormone gene, suggesting gene duplication as a possible cause of the accelerated sequence evolution. Here we determined 21 different GH-like sequences from four species of OWM and hominoids. Combining with published sequences from OWM and hominoids, our analysis demonstrates that multiple gene duplications and several gene conversion events both occurred in the evolutionary history of this gene family in OWM/hominoids. The episode of recent duplications of CSH-like genes in gibbon is accompanied with rapid sequence evolution likely resulting from relaxation of purifying selection. GHN genes in both hominoids and OWM are under strong purifying selection. In contrast, CSH genes in both lineages are probably not. GHV genes in OWM and hominoids evolved at different evolutionary rates and underwent different selective constraints. Our results disclosed the complex history of the primate growth hormone gene family and raised intriguing questions on the consequences of these evolutionary events. © 2005 Elsevier B.V. All rights reserved.

Phylogenetic studies of pantherine cats (Felidae) based on multiple genes, with novel application of nuclear beta-fibrinogen intron 7 to carnivores

The pantherine lineage of the cat family Felidae (order: Carnivora) includes five big cats of genus Panthera and a great many midsized cats known worldwide. Presumably because of their recent and rapid radiation, the evolutionary relationship among pantherines remains ambiguous. We provide an independent assessment of the evolutionary history of pantherine lineage using two complete mitochondrial (mt) genes (ND2 and ND4) and the nuclear beta-fibrmogen intron 7 gene, whose utility in carnivoran phylogeny was first explored. The available four mt (ND5, cytb, 12S, and 16SrRNA) and two nuclear (IRBP and TTR) sequence loci were also combined to reconstruct phylogeny of 14 closely related cat species. Our analyses of combined mt data (six genes; approximate to 3750 bp) and combined mt and nuclear data (nine genes; approximate to 6500 bp) obtained identical tree topologies, which were well-resolved and strongly supported for almost all nodes. Monophyly of Panthera genus in pantherine lineage was confirmed and interspecific affinities within this genus revealed a novel branching pattern, with P. tigris diverging first in Panthera genus, followed by P. onca, P. leo, and last two sister species P. pardus and P. uncia. In addition, close association of Neofelis nebulosa to Panthera, the phylogenetic redefinition of Otocolobus manil within the domestic cat group, and the relatedness of Acinonyx jubatus and Puma concolor were all important findings in the resulting phylogenies. The potential utilities of nine different genes for phylogenetic resolution of closely related pantherine species were also evaluated, with special interest in that of the novel nuclear beta-fibrinogen intron 7. (c) 2005 Elsevier Inc. All rights reserved.

Mitochondrial data support an odd-nosed colobine clade

To obtain a more complete understanding of the evolutionary history of the leaf-eating monkeys we have examined the mitochondrial genome sequence of two African and six Asian colobines. Although taxonomists have proposed grouping the "odd-nosed" colobines

The nuclear matrix of Euglena gracilis (Euglenophyta): A stage of nuclear matrix evolution?

Euglena gracilis cell was extracted sequentially with CSK-Triton buffer, RSB-Magik solution and DNase-As solution. DGD embedment-free electron microscopy showed that in the extracted nucleus there was a residual non-chromatin fibrous network. That it could not be removed by hot trichloroacetic acid further supported the idea that it was a non-histone, non-chromatin fibrous protein network, and should be the internal network of the nuclear matrix. After the sequential extraction, the nuclear membrane was removed, leaving behind a layer of lamina; the chromatin was digested and eluted from the dense chromosomes and residual chromosomal structures that should be chromosomal scaffold were revealed. Western blot analysis with antiserum against rat lamins showed that nuclear lamina of the cell possessed two positive polypeptides, a major one and a minor one, which had molecular masses similar to lamin B and lamin A, respectively. Comparing these data with those of the most primitive eukaryote Archezoa and of higher eukaryotes, it was suggested that the lower unicellular eukaryote E. gracillis already had the nuclear matrix structure, and its nuclear matrix (especially the lamina) might represent a stage of evolutionary history of the nuclear matrix. (C) 2000 Editions scientifiques et medicales Elsevier SAS.

The complete mitochondrial DNA sequence and the phylogenetic position of Achalinus meiguensis (Reptilia: Squamata)

The mitochondrial genome complete sequence of Achalinus meiguensis was reported for the first time in the present study. The complete mitochondrial genome of A. meiguensis is 17239 bp in length and contains 13 protein-coding genes, 22 tRNA, 2 rRNA, and 2 non-coding regions (Control regions). On the basis of comparison with the other complete mitochondrial sequences reported, we explored the characteristic of structure and evolution. For example, duplication control regions independently occurred in the evolutionary history of reptiles; the pseudo-tRNA of snakes occurred in the Caenophidia; snake is shorter than other vertebrates in the length of tRNA because of the truncations of T psi C arm (less than 5 bp) and "DHU" arm. The phylogenic analysis by MP and BI analysis showed that the phylogenetic position of A. meiguensis was placed in Caenophidia as a sister group to other advanced snakes with the exclusion of Acrochordus granulatus which was rooted in the Caenophidia. Therefore we suggested that the subfamily Xenodermatinae, which contains A. meiguensis, should be raised to a family rank or higher rank. At the same time, based on the phylogenic statistic test, the tree of Bayesian was used for estimating the divergence time. The results showed that the divergence time between Henophidia and Caenophidia was 109.50 Mya; 106.18 Mya for divergence between Acrochordus granulatus and the other snakes of the Caenophidia; the divergence time of A. meiguensis was 103 Mya, and Viperidae diverged from the unilateral of Elapidae and Colubridae was 96.06 Mya.

Bioinformatic identification of genes encoding C1q-domain-containing proteins in zebrafish

C1q is the first subcomponent of classical pathway in the complement system and a major link between innate and acquired immunities. The globular (gC1q) domain similar with C1q was also found in many non-complement C1q-domain-containing (C1qDC) proteins which have similar crystal structure to that of the multifunctional tumor necrosis factor (TNF) ligand family, and also have diverse functions. In this study, we identified a total of 52 independent gene sequences encoding C1q-domain-containing proteins through comprehensive searches of zebrafish genome, cDNA and EST databases. In comparison to 31 orthologous genes in human and different numbers in other species, a significant selective pressure was suggested during vertebrate evolution. Domain organization of C1q-domain-containing (C1qDC) proteins mainly includes a leading signal peptide, a collagen-like region of variable length, and a C-terminal C1q domain. There are 11 highly conserved residues within the C1q domain, among which 2 are invariant within the zebrafish gene set. A more extensive database searches also revealed homologous C1qDC proteins in other vertebrates, invertebrates and even bacterium, but no homologous sequences for encoding C1qDC proteins were found in many species that have a more recent evolutionary history with zebrafish. Therefore, further studies on C1q-domain-containing genes among different species will help us understand evolutionary mechanism of innate and acquired immunities.

Utility of ITS1-5.8S-ITS2 sequences for species discrimination and phylogenetic inference of two closely related bucephalid digeneans (Digenea : Bucephalidae): Dollfustrema vaneyi and Dollfustrema hefeiensis

The complete internal transcribed spacer 1 (ITS1), 5.8S ribosomal DNA, and ITS2 region of the ribosomal DNA from 60 specimens belonging to two closely related bucephalid digeneans (Dollfustrema vaneyi and Dollfustrema hefeiensis) from different localities, hosts, and microhabitat sites were cloned to examine the level of sequence variation and the taxonomic levels to show utility in species identification and phylogeny estimation. Our data show that these molecular markers can help to discriminate the two species, which are morphologically very close and difficult to separate by classical methods. We found 21 haplotypes defined by 44 polymorphic positions in 38 individuals of D. vaneyi, and 16 haplotypes defined by 43 polymorphic positions in 22 individuals of D. hefeiensis. There is no shared haplotypes between the two species. Haplotype rather than nucleotide diversity is similar between the two species. Phylogenetic analyses reveal two robustly supported clades, one corresponding to D. vaneyi and the other corresponding to D. hefeiensis. However, the population structures between the two species seem to be incongruent and show no geographic and host-specific structure among them, further indicating that the two species may have had a more complex evolutionary history than expected.

Sox genes in grass carp (Ctenopharyngodon idella) with their implications for genome duplication and evolution

The Sox gene family is found in a broad range of animal taxa and encodes important gene regulatory proteins involved in a variety of developmental processes. We have obtained clones representing the HMG boxes of twelve Sox genes from grass carp (Ctenopharyngodon idella), one of the four major domestic carps in China. The cloned Sox genes belong to group B1, B2 and C. Our analyses show that whereas the human genome contains a single copy of Sox4, Sox11 and Sox14, each of these genes has two co-orthologs in grass carp, and the duplication of Sox4 and Sox11 occurred before the divergence of grass carp and zebrafish, which support the "fish-specific whole-genome duplication" theory. An estimation for the origin of grass carp based on the molecular clock using Sox1, Sox3 and Sox11 genes as markers indicates that grass carp (subfamily Leuciscinae) and zebrafish (subfamily Danioninae) diverged approximately 60 million years ago. The potential uses of Sox genes as markers in revealing the evolutionary history of grass carp are discussed.